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ON THE ASYMPTOTIC DISTRIBUTION OF BLOCK-MODIFIED 

RANDOM MATRICES 

OCTAVIO ARIZMENDI, ION NECHITA, AND CARLOS VARGAS 


Abstract. We study random matrices acting on tensor product spaces which have been trans¬ 
formed by a linear block operation. Using operator-valued free probability theory, under some mild 
assumptions on the linear map acting on the blocks, we compute the asymptotic eigenvalue distribu¬ 
tion of the modified matrices in terms of the initial asymptotic distribution. Moreover, using recent 
results on operator-valued subordination, we present an algorithm that computes, numerically but 
in full generality, the limiting eigenvalue distribution of the modified matrices. Our analytical re¬ 
sults cover many cases of interest in quantum information theory: we unify some known results and 
we obtain new distributions and various generalizations. 
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1. Introduction 

We continue the work of Banica and second named author [BN12b] on finding the asymptotic 
distribution (as d —>• oo) of dm x dm random matrices which have been block-wise modified by 
some fixed self-adjoint linear transformation (p : M m ( C) —»• M n { C). More concretely, we consider a 
self-adjoint dm x dm unitarily invariant random matrix Xd and a linear map p : M m ( C) —>• M n (C). 
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Our goal is to understand the asymptotic eigenvalue distribution of the block-wise modified random 
matrix 

X% := [idd ® <p](X d ) E M d ( C) ® M n ( C). 

The motivation for this line of research comes from quantum information theory. The problem 
of deciding whether a quantum state p E M n (C) <g> M m (C) is entangled has been shown to be 
equivalent to the fact that, for all positivity preserving maps ip : M m ( C) —> M n ( C), the modified 
state p = [id <8> ip] (p) is positive semi-definite [HHH96]. It is thus important to understand how the 
spectrum of generic matrices behaves under linear block-transformations. Aubrun [Aubl2] studied 
the positivity of partially transposed random quantum states: this corresponds to the transposition 
map acting on the blocks of a Wishart element. He computed the range of parameters of the 
Wishart distribution for which the partial transposition is still positive. Similar computations, in 
different asymptotic regimes, or for different maps acting on the blocks, have been performed in 
[BN12a, FSl3, JLN14, JLN15]. In this paper, we unify these results by showing that operator¬ 
valued free probability theory provides the right framework to study such questions. 

Our new approach relies on the theory of operator-valued free probability, developed by Voiculescu 
[Voi85, Voi95] and later by Speicher [Spe98]. We find a general numerical solution using the tools 
from [BSTV14] which allow to approximate operator-valued free multiplicative convolutions. For 
certain cases, we are able to get explicit formulas for the distributions by using the operator-valued 
S-transform [Dyk06] over some suitable commutative algebras. 

For the case m = n, the main observation is that we may write 

m 

X d= E uZ(I d ®E i:j )X(I d ®E kl ), 

i,j,k,l =1 

where E tJ E M m (C) are the matrix units. The joint non-conmiutative distribution (see Section 2 
for the main definitions) of the matrices (I d (g> -Fy )f) =1 (with respect to the state Tdm :=Eo g^Tr) 
does not change with d. Thus, if the random matrix X is unitarily invariant (for example, if X 
is a Wishart, Wigner, or randomly rotated matrix), we are allowed to use Voiculescu’s asymptotic 
freeness results [Voi91] to compute the asymptotic joint moments of the non-commutative random 
variables X , ( I d (g> A,j)™ =1 as d —> oo. These asymptotic joint moments are given by replacing 
(Id <8> Eji)™ = i and X by an abstract collection of operators (ey)ij< m , x in some non-commutative 
probability space (A,t). 

The limiting distribution of X ^ is the same as the distribution of the element 

m 

x v ■- a^e ijX e k i. 

i,j,k,l =1 

Thus, the problem is reduced to study the distribution of such elements for which the machinery of 
operator-valued free probability is available. Note that the general case m / n is solved similarly 
with the aid of rectangular spaces, using the theory developed by Benaych-Georges [BG07, BG09a, 
BG09b] . 

The paper is organized as follows. We start with a short overview of operator-valued free prob¬ 
ability, the framework in which we state all our results. Section 3 contains the numerical solution 
to the problem, in full generality. We then move to the theoretical study: in Section 4 we prove 
a freeness result, under some assumptions on the eigenvectors of the Choi matrix C of the linear 
map ip acting on the blocks. The exact distribution of the modified matrix is computed in Section 
5, under more stringent assumptions on C. Sections 6 and 7 contain some important examples, 
both of distributions and of maps ip ; most of the examples are of interest in Quantum Information 
Theory. 
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2. Operator-valued free probability 

Free probability was introduced by Voiculescu [Voi85] in order to tackle some problems in op¬ 
erator algebras. Later, the theory detached from its operator-algebraic origins, and grew into a 
branch of non-commutative probability, where the concept of freeness replaces the notion of in¬ 
dependence from classical probability. The main object of study are (non-commutative) random 
variables a\,... ,dk in a (V-algebra A and their joint distribution with respect to a given state 
(i.e. a unital positive linear functional r : A —> C). The pair (A,t) is called a non-commutative 
probability space. 

By the joint distribution of an ordered tuple a = (a \,..., a*;) of random variables in (A, t), we 
mean the collection of maps <L(a) = Ur>o > where 

<f>“ := <F r = {(ii,..., i r ) r(a il ... a ir ) : i u ..., i r < k} 

are the r-th order mixed moments of (ai,..., aif). 

The state t should be thought as playing the role of the expectation functional in classical prob¬ 
ability. Motivated by operator-algebraic constructions, Voiculescu defined freeness in [Voi85], as a 
new, non-commutative notion of independence. In [Voi91], free random variables where constructed 
as limits of certain large random matrices, establishing a connection between free probability theory 
and random matrix theory. 

The idea of operator-valued free probability ( [Voi95] ) is to generalize free probability, by replacing 
the scalar valued linear form t : A —> C by a conditional expectation E : A —> B onto a larger 
sub-algebra C C B C A. This leads to a broader definition of freeness, which occurs in more general 
situations of random matrix theory, as observed already by Shlyakthenko in [Shl96] . 

Many aspects of the theory of (scalar-valued) free probability can be lifted to the operator¬ 
valued level. The combinatorics of operator-valued free probability (see [Spe98]) remains the same 
provided that the nesting structure of non-crossing partitions is respected while operating. This 
will make cumulants particularly handy when finding good candidates for the smallest sub-algebra 
B C A over which two given random variables are free. 

The question of finding explicit formulas for operator-valued distributions is closely related to 
the possibility of finding realizations of the distribution in terms of operators which are free over a 
commutative algebra. We use the criteria in [NSS02] to give sufficient conditions (in terms of the 
Choi matrix of the map ip) for such a realization to exist. 

For the cases where B is non-commutative it is extremely complicated to obtain exact distribu¬ 
tions. On the other hand, numerical algorithms which rely on subordination, such as the ones in 
[BB07] for the computation of the additive and multiplicative free convolutions admit very effective 
generalizations to the operator-valued level (see [BSTV14, BMS13]). We will use these to obtain a 
general numerical solution to our problem. 

2.1. Basic definitions. We gather here some basic definitions that will be needed in what follows, 
and we establish some notation. 

Definition 2.1. (1) Let A be a unital *-algebra and let C C B C A be a *-subalgebra. A B- 
probability space is a pair (A, E), where E : A B is a conditional expectation, that is, a linear 
map satisfying: 

E (bob') = bE(a)b', V6, b' eB,a£ A 
E (1) = 1. 

(2) Let (A, E) be a B-probability space and let a := a — E(a)l^ for any a € A. The *-subalgebras 
B C A \,..., Ah C A are B-free (or free over B, or free with amalgamation over B) (with respect to 

K)iff 


E(aia 2 • • • a r ) = 0, 


(1) 
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for all r > 1 and all tuples ai,...,a r G A such that ai G AjU) with j( 1) ^ j{2) ^ ^ j(r). 

(3) Subsets S\,..., Sk C A are B-free if so are the *-subalgebras (S i, B ),..., (S^, i3). 

Similar to independence, freeness allows to compute mixed moments free random variables in 
terms of their individual moments. 

The most basic examples of non-commutative probability spaces are the space of deterministic 
matrices (M n ( C), -TV) and the space of classical random variables (A, IE). These serve as building 
blocks for more general (operator-valued) non-commutative probability spaces. 

Example 2.2 (Rectangular probability spaces). Let A be a unital C*-algebra and let r : A -A C 
be a state. Let pi,. ■ ■ ,Pk be pairwise orthogonal projections with 1 = p\ + • • • + Pk and T(pf) > 0 
for all i. There exist a unique conditional expectation E : A -A (p±,... ,pk) compatible with r in the 
sense that r o E = r. The conditional expectation is explicitly given by 

E(a) = y, T(p i )~ 1 T(p i api). 
i<k 

Definition 2.3. Let (An, tjv), N > 1, and(A,r ) be C -probability spaces and let (a[ N \ ... , aj^) G 
A%, (ai,..., a*) G A k be such that 

Jim r N ((a\* ] ) ■ ■ • (a[f } )) = r(a h • • • a ir ), 

iv—>-oo ' 

for all r > 1, 1 < ii,... ,i m < k Then we say that {a[ N \ ..., aj^) converges in distribution to 
(ai,..., afc) and we write (a[ N \ ..., a^) -A (ai,..., a*,). 

Remark 2.4. (1) The joint distribution (with respect to T m d) of the blown-up matrix units ( I c i ® 
Eij)i,j<m Is independent of d. We may consider an abstract “limiting” non-commutative probability 
space {A,t) such that M m ( C) C A and t\m/q = ^Tr. If ( eij)i,j< m denote the matrix units in 
M m (C) C A we have that 

(En <8> Id, E 12 ® Id,..., E mm (8> Id) -A (en, ei 2 , • • ■ j &mm), 
with joint distribution determined by the usual rules: 

m 

f'ij&kl = djkeu, r(eij) = n dij, ^ ^ £a = 1) = &ji- 

i =1 

(2) If Xd is an mdxmd Wigner, Wishart or randomly rotated self-adjoint matrix with lim^oo Td n (X(j) 
r(x k ) for some fixed operator x in a non-commutative probability space (A,t), then by [Voi91], we 
have that 

(X, En ® Id, E \2 <8 Id, ■ ■ ■, E mm <8 Id) t (x, eu , ei 2 , • • •, e mm ), 
where x and (eij)ij< m are free. This determines completely the joint distribution of these variables. 

(3) Let us now replace m by m + n in the first part of this remark so that our limiting non- 
commutative probability space (A,t) contains now a copy o/(M( m+n J C), ^Jr^Tr). If Xd is still an 
md x md random matrix, as in part (2) (completed with zeros to an d(m + n) x d(m + n) matrix), 
then we have again, that 

(Xd, En <8 Id, E\2 <8 Id, ■ ■ ■ , ^'(m+n)(?n+n) *8 Id) t (x, en, ei2, • • ■ , C( m +n)(m+n))• 

It is quite easy to see that x and ((e?//)*,i<m+n) are no longer C -free. Indeed, the joint distribution 
is aware of the orthogonality relations between Xd and the different matrix units, and we have, for 
example, that 

0 t(x e( m _|_i)( m+ i)) A t(x )T(e( m _|_i)( m _|_i)), 
and hence freeness is broken unless x = 0. 
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As noticed in [BG09b] ; this problem gets fixed by considering distributions with values on the 
commutative algebra generated by the projections 

P<i,m ■ ^ ' Id, ® I'll ■ Pd,n ■ /* ' Id ® I'l i ■ 

i<m m<i<m-\-n 

which in turn converge in distribution to the orthogonal projections 

m ra+n 

Pm • = ^ ^ &a and p n '.= ^ ^ eH) 

i=l i=m -\-1 

with non zero traces r(p m ) = m/(m + n), r(p n ) = n/(m + n). Then x and (( e «j)i,j<m+n ) are 

( Pm, Pn)-free. 

Example 2.5 (Matrix-valued probability spaces). Let A be a unital C*-algebra and let t : A —> C 
be a state. Consider the algebra M n (A) = M n ( C) < 8 > A of n x n matrices with entries in A. The 
maps 

E 3 • (aij)ij 1 y (t (o,{j ))ij G M n (C), 

T 2 ■ (a tJ )jj i >■ ( bjj t{ djj j )jj G D n (C), 

and 

n 1 

Ei . (o,jj fj i ^ ^ ' T((ijj)/ ri G C • I n 

i =l 

are respectively, conditional expectations onto the algebras M n ( C) D D n (C) D C • I n of constant 
matrices, diagonal matrices and multiples of the identity. 

If A\,..., A}, are free in (M, r), then the algebras M n (Ai ),..., M n (Ak) of matrices with entries 
in Ai ,..., Ak respectively are in general not free over C (with respect to Ei). They are, however 
M n (C)-free (with respect to E 3 ). Below is a slightly more general assertion of this simple but 
fundamental result. 

Proposition 2.6. Let (A, E) be a B-probability space, and consider the M n (B)-valued probability 
space (M n ( C) <S> A, id® E). If A \,..., A^ C A are B-free, then (M n (C) < 8 > Ai),..., ( M n (C ) (g> Ajf) C 
(M n (C)<g>A) are ( M n (B))-free. 

Proof. Let , a d 1 ) G M n (C ) (g> A be such that ad G M n ( C) (g> Aj^ with j( 1) / j(2) A ■ ■ ■ A 

j(m). Observe that 

dL = ad - (id®E)(a«) = ((ad) - E(a«)) rs < n = (a#)„< n . 

Hence 

n _ _ _ 

(«®E)((^))...(S^)))= Y. •••(“£,<„)))«<„ =0. (2) 


2.2. Combinatorics and cumulants. A partition of a set is a decomposition into disjoint subsets, 
called blocks. The set of partitions of the set [n] := {1,... ,n} is denoted by V(n). Any partition 
defines an equivalence relation on [n] and vice versa: given n G V(n), i ~ 7r j holds if and only if 
there is a block V G vr such that i,j G V. 

A partition ir G V(n ) is non-crossing if there is no quadruple of elements 1 <i<j<k<l<n 
such that i k, j l and i A™ j • The non-crossing partitions of order n form a sub-poset of 
V ( n ) which we denote by J\fC (n). 






6 


O. ARIZMENDI, I. NECHITA, AND C. VARGAS 


For n G N, a C-multi-linear map / : A n —>■ £> is called B-balanced if it satisfies the £>-bilinearity 
conditions, that for all b,b' E B, a\,... ,a n E A, and for all r = 1 ,...,n — 1 

/ (bai,... ,a n b') = bf (ai,..., a n ) b' 
f (o>i, • • •, a r b, CL r -\~ i, • • •, a> n ) — / (oi,..., a r , ba r -fi ..., a n ) 

A collection of £>-balanced maps (fn) w& jyc is said to be multiplicative with respect to the lattice 
of non-crossing partitions if, for every n G A fC, f n is computed using the block structure of it in 
the following way: 

• If 7 T = l n G MC (n), we just write f n := f n . 

• If l n f n = {V \,..., 14} G MC ( n ), then by a known characterization of MC, there exists 
a block V r = {s + l,...,s + /} containing consecutive elements. For any such a block we 
must have 


fn (®1 1 • • • ) Q>n) — fn\V r ( a l> ■ • ■ > Q>sfl ( a s+l j ■ ■ ■ > a s+i) i &S+/+1 j • • • > a n) j 
where 7 r\I 4 G AAC (n — l) is the partition obtained from removing the block V r . 

We observe that a multiplicative family (fn) ne j\rc entirely determined by (/n) ngN . On the other 
hand, every collection (/ n ) ngN of £>-balanced maps can be extended uniquely to a multiplicative 
family (fn) neJ vc- 

The operator-valued free cumulants (-Rf) 7re _^ c are indirectly and inductively defined as the 
unique multiplicative family of £>-balanced maps satisfying the (operator-valued) moment-cumulant 
formulas 

IE (cii ... u n ) — ^ ) R n (oi, • • •, a n ) 


n£jVC(n) 

By the cumulants of a tuple a\, ..., G A, we mean the collection of all cumulant maps 


R 


JS]CL\ 

^1 5 ‘“An 


B n ~ l -a B, 

(bi,...,b n - 1 ) i-A R% (ai 1 ,biai 2 ,...,bi n _ 1 a in ) 


for n G N, 1 < ii, ..., i n < k. 

A cumulant map Rj’ ai ’y ,ak is mixed if there exists r < n such that i r A i r + 1 - The main feature 
of the operator-valued cumulants is that they characterize freeness with amalgamation: 


Theorem 2.7 ([Spe98]). The random variables a\,... ,a n are B-free iff all their mixed cumulants 
vanish. 


We recall a useful result from [NSS02] which gives conditions for (operator-valued) cumulants to 
be restrictions of cumulants with respect to a larger algebra. 


Proposition 2.8. Let lEVcBcAbe algebras such that (A, IF) and (B, E) are respectively 
B-valued and V-valued probability spaces (and therefore (A, IE o F) is a V-valued probability space) 
and let ai,... ,ak E A. 

Assume that the B-cumulants of oi,..., a*, G A satisfy 

i^::x ah (di,...,d*-i)GP, 


for all n E N, 1 < ii, ■.., i n < k, d \,..., d n -\ G V. 

Then the V-cumulants of aq,..., a*, are exactly the restrictions of the B-cumulants of aq,..., a* 
namely 


/or fllinGN, 1 < *i,... An < k, d \,..., d n -i £ T>. 




(dl,. . . , d n —1 ) , 


Corollary 2.9. Let 13 C Ai, A .2 C A 6e B-free and let T> C M/v(C)®£>. Assume that, individually, 
the Mjv(C)<8> B-valued moments (or, equivalently, the M^(C) <g>B-cumulants) of both x G Mjv( C)<8> 
Aq and y E Mjv(C) <8 > A 2 , rcLen restricted to arguments in V, remain in V. Then x, y are D-free. 
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Proof. By Proposition 2.6 x,y are Mjv(C) ® B free. The Mjv(C) <g> B-mixed moments on 
restricted to arguments d\,..., d n -\ G P can be expressed through the Mjy(C) ® £> moment- 
cumulant formula 

E (aidi... d n -\a n ) = ^ Pf’ v(C) ® B (ai,..., dn-PP), 

n£j\[C(n) 

as a sum of (nested) products of individual Mjv(C) <8> £> cumulants on ar and y. For each partition 
7 r, we pick an interval block. The cumulant corresponding to this block remains in P (in particular 
it is zero if the cumulant is mixed). By proceeding inductively removing interval blocks, we obtain 
that each sumand is in V, and hence, so is any restricted mixed moment E (oidi... d n -\a n ). 

By Moebius inversion (or induction), any restricted mixed cumulant is also in V and we are 
allowed to use Theorem 2.8. Hence the mixed T> cumulants of x,y coincide with the restrictions of 
the Mjv(C) <S)B cumulants, which vanish due to (C) ® B-freeness. Thus the P-mixed cumulants 
vanish as well, proving the assertion. □ 

In view of Remark 2.4 (3), if A\ := ((e lJ ),; J <m+n) and A 2 := (x) are ( Pm,Pn ) free, then (by 
Proposition 2.6) M r { C) <8> A\ and M r ( C) (g> A -2 are M r { C) ® (p m ,Pn}- free. Our goal is to replace 
M r ( C) ® (pm,Pn )-freeness by £>-freeness for some commutative subalgebra B C M r { C) ® (p m .■ Pn ), 
using Corollary 2.9. 

2.3. Transforms. Like in the scalar-valued case, there are analytical tools to compute operator¬ 
valued free convolutions, which are based on the B-valued Cauchy transform 

G'f(6) = E((6-x)- 1 ), 

which maps the operatorial upper half-plane 

H + (£>) := {b G B : Be > 0 such that — i(b — b*) > e • 1} 

into the lower half-plane H _ (B) = —H + (B). In the usual settings coming from random matrix 
models (as we have seen above), our probability space A may have several operator-valued structures 
Ej : A -A B{ simultaneously, with C = B\ C £>2 C ••• C Bk, and Ej o Ej+i = Ej. We are 
usually interested ultimately in the scalar-valued distribution, which can be obtained (via Stieltjes- 
inversion) from the Cauchy transform. The later can be obtained from any ’’upper” £>,;-valued 
Cauchy transform, as we have that, for all b G Bj 

E i (Gx i+1 (b))=E i oE i+1 ({b-x)- 1 )=E i ({b- X )~ 1 ) = G^(b). 

In terms of moments, the operator-valued Cauchy transform is given by 

G*(b) = E((6 - x)- 1 ) = £ E (b- l {xb~ l ) n ). 

n> 0 

The operator-valued P-transform is defined by 

( 6 ) = bx ' ■ ■ ■ >• 
n>l 

The vanishing of mixed cumulants for free variables implies the additivity of the cumulants, and 
thus also the additivity of the P-transforms ([Voi95]): If a± and 02 are B-free then we have for 
6 G B that Pf 1+Q2 (6) = Pf x (6) + P a2 (6). 

These transforms satisfy the functional equation 

Gf(6) = (Pf (Gf(6))-6) _1 . (3) 

The 5-transform was dehned by Voiculescu [Voi87] to compute multiplicative free convolutions. 
The most general operator valued generalization of the S transform was defined in [Dyk06]. 
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If M a (b) = J2k> i IE((a. 6 ) fc ) is the moment generating series of a, it is possible to define a com¬ 
positional inverse Afi 1 ^( 6 ) in a suitable domain using the analytic inverse function theorems on 
Banach algebras. The 5-transform is then defined as S a (b) = (1 + b)b~ 1 Ma b ). 

The 5-transform factorizes the operator-valued multiplicative convolution in the following sense: 
If x, y are 5-free, then 

S xy (b) = S y (b)S x (S y (b)- l bS y (b)). 

In this work we will be mainly interested in the case where the algebra 5 is commutative, the 
right hand can be then simplified to S x (b)S y (b). 

For our general numerical solution, where B is not commutative, the factorization property of 
the 5-transfornr is implicitly being used in the subordination formulas (see [BSTV14]). 

In the commutative case, we will use the following equation relating the 5 and the ^-transforms: 
blZ(b) +S(blZ(b)) = b (see, for example, Lecture 18 in [NS06]). 

A drawback of the operator-valued setting is that, unless we ask B to be commutative, one 
can hardly compute explicit distributions: although we have 5-valued generalizations of the 1Z 
and 5-transfornrs, the task of explicitly inverting these operator-valued analytic maps is nearly 
impossible for any non-trivial situation (even for finite dimensional, relatively simple sub-algebras, 
like B = M 2 (C)). 

A very powerful (numerical) method to obtain 5-valued free convolutions is based on the analytic 
subordination phenomena observed by Biane ([Bia98], see also [Voi02]). In particular, the approach 
of [BB07] to obtain the subordination functions by iterating analytic maps can be very efficiently 
performed in the 5-valued context. 

In [BSTV14] the fixed point algorithm from [BB07] for computing free multiplicative convolutions 
was generalized to the operator-valued level. Let us first recall the eta transform (or Boolean 
cumulant transform) 


Vx(b) = 1 - b(E((b 1 -x ) x )) 1 . 

Theorem 2.10 ([BSTV14]). Let x > 0, y = y* £ A be two random variables with invertible 
expectations, free over B. There exists a Frechet holomorphic map oj 2 : {b £ 5: Zs(bx) > 0} —> 
H + (S), such that 

( 1 ) g y (uj 2 (b)) = r/ X y(b), Zs(bx) > 0 ; 

( 2 ) ui 2 {b) and b~ l ui 2 {b) are analytic around zero; 

(3) For any b £ 5 so that Zs(bx) > 0, the map g H + (B) -A H + (5), gb(w) = bh x {h y (w)b), 
where 

h x {b) = 6 _1 - E [( 6 -1 - x)^ 1 ] 1 - (4) 

is well-defined, analytic and for any fixed w G H + (5), 

oj 2 (b) = lim gl n (w), 

n—>oo 

in the weak operator topology. 

Moreover, if one defines coi(b) := h y (co 2 (b))b, then 

Vxy(b) = w 2 (%z(wi(&))w 2 (fr) - \ Zs{bx) > 0. 

The invertibility condition on E(y) can be dropped if we restrict to finite dimensional algebras. 

Proposition 2.11 ([BSTV14]). Let B be finite-dimensional and x > 0, y = y* be free over B. 
There exists a function g: {b G 5: Zs(bx) > 0} x H + (5) -A H + (5) so that 

(1) ui 2 {b) := limn-xx, gl n (w) exists, does not depend on w G H + (5) J and is analytic on {b G 
B:Q{bx)> 0}; 
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( 2 ) 

r] y (u2(b )) = r} xy (b), be {be B: %(bx) > 0 }. 

From a numerical perspective, switching between h x to G x is a simple operation. By Proposition 
2.11 one only needs the individual B -valued Cauchy transforms of x,y (or good approximations of 
these) to obtain the B-valued Cauchy transform of xy, and hence, its probability distribution. 

In a scalar-valued non-commutative probability space (A, t), we have the integral representation 
of the Cauchy transform: 

G x (z) = t((z - xy 1 ) = [ (z -t)~ l dy x (t). 

Jr 

Analogously, for linear, self-adjoint elements c® x in a M n (C)-valued probability space (M n (C) < 8 > 
(A), id m < 8 ) r), we have: 

Gc® x {b) = (id m <g> t)((6 - c<g> x)^ 1 ) = / (b - c®ty l dn x (t). (5) 

Jr 

The previous integrals can be approximated, for example, by using matrix-valued Riemann sums. 

The case of deterministic matrices is simpler. If we assume that M n (C) C A and consider 
x = x* e M m { C) <S> M n ( C). Then G x (b) = G x (b ® I n ) is just the partial trace of the resolvent 

G x (b) = (id m <g> -Tr)(( 6 <g>/ n -x)^ 1 ). ( 6 ) 

n 

2.4. Unitarily invariant random matrices: definition and examples. In this work, our main 
focus is the effect of a block-linear transformation on the (asymptotical) spectrum of a random 
matrix. The models of random matrices we consider are as follows. 


Definition 2.12. A self-adjoint random matrix X e M n (C) is called unitarily invariant if, for all 
unitary operators U eU n , the random matrices X and JJXU* have the same distribution. 


Let us consider two special cases of the definition above: the GUE ensemble and the Wishart 
ensemble (see [AGZ10] for the general theory of random matrices). 

Both GUE and Wishart ensembles are built out of i.i.d. complex Gaussian random variables. For 
any integer size parameter n, consider the self-adjoint random matrix G n having upper diagonal 
entries given by 


G n {i,j) 


’An / \fn, 

( Zij + 4)/V2n, 


if i = j 

if i < j, 


where Zij and Zij are i.i.d. standard complex Gaussian random variables. The classical Wigner 
theorem [AGZ10, Theorem 2.2.1] states that the non-commutative random variables G n converge 
in distribution to the standard semicircle distribution 


yj 4- - 

/Uc; 0,1 = - Y~ -1 [— 2 , 2 ] (x)dx. 

More generally, the semicircle distribution with average a e M and variance a 2 > 0 is defined by 

_ y / 4u 2 - (x-a) 2 ^ 

/^sc;a,cr — 27T(J^ l[a— 2cr,a+2a] \%)d>X. 

The Wishart ensemble is defined starting from a rectangular random matrix X n e M nX k{ C), 
having i.i.d. standard complex Gaussian entries. A Wishart matrix of parameters (n, k) is then 
W n = X n Xf. Consider now, for a sequence of parameters k n such taht k n /n —> c > 0, a sequence 
of random Wishart matrices W n of parameters ( n,k n ). The limiting eigenvalue distribution of 
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the sequence W n , as n —> oo, is given by the Marcenko-Pastur (or free Poisson) distribution of 
parameter c (see [AGZ10, Exercise 2.1.18]) 


/ x c- \/4c — (x — 1 — c) 2 
7r c = max(l — c, 0)5 q H---1 


27 TX 


l [1+c-2V5,1+c+2V5] 


dx. 


The family of free Poisson measures 7r c are characterized, in terms of free probability, by the fact 
that their free cumulants are all equal to c. This family admits the following extension. For 
a compactly supported probability measure i , define the compound free Poisson distribution of 
parameter /a as the unique probability measure having free cumulants given by 


Vn > 1, 



x n dp(x). 


Finally, let us introduce the model of unitarily invariant random matrices we are interested 
in, and of which both the GUE and the Wishart ensembles are examples of. For a compactly 
supported probability measure on the real line fi, we associate a sequence of self-adjoint random 
matrices X n £ M n ( C) defined as 

X n = U n D n U*, 

where D n is a sequence of deterministic, diagonal matrices converging in distribution to /r, and 
XJ n £ U n is a sequence of Haar-distributed random unitary matrices. In this paper, we are going 
to suppose that the space on which X n acts has a tensor product structure, and we are going to 
study the eigenvalue distribution of the modified version of X n , obtained by acting with a fixed 
linear map p on the blocks of X n . 

In view of Remark 2.4, in order to understand the asymptotic distribution of a block-modified 
unitarily invariant random matrix ensemble, it seems important to study first the asymptotic joint 
distribution of its blocks (xij), where Xij = euxej\. By Remark 2.4, this joint distribution depends 
only on m and p. 

In Lecture 20 of [NS06] it was observed that, with respect to the compressed state Tn(a) = 
mr(enaen), the joint distribution of ( x ^) is completely characterized by their mixed cumulants, 
which satisfy the R-cyclic property: 


W 1 1 


c hji i x i23i ; 


, X■ 


't'kjk ) 


m (fc ^r k (x,x,...,x), if ji = h, 32 = k, ■ ■ ■, jk = k, 
0 , otherwise 


( 7 ) 


where r n (x,x,... ,x) are the free cumulants of x (or n) in the original space. This implies that, 
more generally, for any 7r £ J\fC{k ), 

• • ,Xi k ik) = m~ {k ~ M) r^(x,x, ...,x), 


whenever the cyclic condition j Vl = i V2 ,... ,j v ^ = i Vl holds for each block V = {ui,..., U|y|} £ ir, 
and zero otherwise. 

The distribution of a block-modified ensemble will then be that of the element 


x* := ^2 °ij ® Xij £ M n ( C) <g) -An 

i,j<m 


with respect to the functional := tr n <g> m, where Cij = p>(e{j) are the blocks of the Choi matrix 
of p> : M m { C) —> M n { C) (see Section 4). 

For example, the distribution of each diagonal sumand cu <S> xn = c* u ® x is given by the 
classical multiplicative convolution of the discrete measure n Cii with the distribution p Xll of the 
free compression of x. Furthermore, the R-cyclic condition (7) implies that the mixed cumulants 
of diagonal elements vanish. Hence in,... ,x mm 6 An are identically distributed and free among 
themselves. 
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For a general situation we get that the moments of are given by 


tp((xV) k ) = ® X ilh ) 

■ ■ ■ ( <H k j k ® x ikjk )) 

(8) 




= t r n(ci 1 j 1 ■ ■ ■ Ci k j k 

) T ll( x iljl ■ ■ ■ Xi k j k ) 

(9) 


h,ji,—,ikdk<n 

= tr n (ci 1 j 1 ■ ■ ■ Ci k j k ) R n 11 (xi 1 j 1 ■ ■ ■ Xi k j k ). (10) 

h Tr£j\fC(n) 

Hence the joint distribution of the Choi blocks and the blocks of unitarily invariant random 
matrices interact in a non-trivial way. One of the main objectives in this paper is to find conditions 
on the Choi matrices so that we are able to compute the distribution of x v explicitly. 


3. General numerical solution 

Now we show that the free multiplicative convolution can be used to compute the distribution of 
random matrices which have been deformed block-wise by a fixed self-adjoint linear transformation. 
For simplicity and clarity of exposition, we only solve the case n = m (from our treatment of the 
case to 7 ^ n in the subsequent Sections, it should be clear how to adapt our algorithm to cover the 
general case). 

We have already pointed out that the limiting distribution of Xf is the same as the distribution 
of 

m 

xv -■= a ki e u xe kh 

i,j,k,l= 1 

where are matrix units, free from x, and aV, are fixed constants. The self-adjointness 

condition on (p is equivalent to the condition a' kl = d 1 ^ , for all i,j,k,l < to. 

If <p is self-adjoint, it will be useful to express it as a difference of positive maps. For this we 
consider the lexicographic order of two-letter words in an m-letter alphabet, and write, for each 
(i,j) < (, l,k ), 

re t9 eijxe k i + re~‘ e ei k xeji = r l/2 (e l9/2 eij + e~ l6/2 ei k )xr 1/2 (e~ l9/2 eji + e l9/2 e k i ) 

reijxeji tci^xc^i 

= (fkl) x (fkl)* ~ re ij xe ji ~ re lk xe kl 

so that we get 

X<P = J2 ( fklXfkl)* + J2 Pijeijxeji 

l<i,j,k,l<m l<i,j<m 

( i,j)<(l,k ) 

= E ifSwiir+ E /«(-i) eW W« 

(i,j)<(l,k) 

J2 (fki) £ ki x (fkiY > 

l<i,j,k,l<m 

(i,j)<(l,k) 

where eY = (—l) £ (*>i) and ^ = \ for (i,j) / (l, k). 

From the elements /(,], (i,j) < ( l,k ) we build a vector / = (,/'n , f \2 , fmm ) °f size N := 

to 2 (to 2 + l)/2. We consider also the diagonal matrix E = diag(£}},£}),...,£“) and we let 
x := E <g> x. 
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We see that fxf* = x^, so the desired distribution is the same (modulo a Dirac mass at zero of 
weight 1 — 1/N) as the distribution of f* fx in the C *-probability space (M/v(C) 0 A, -^Tttv 0 r). 
Moreover, since w and each of the e tJ are free, by Proposition 2.6, the matrices f*f and x are free 
with amalgamation over B = Mjy( C)0l (with respect to the conditional expectation E := idj\r0r). 

By Proposition 2.11, we can obtain the Mtv(C)- valued Cauchy transform of fxf* numerically, 
provided that we can compute the M 7 v(C)-valued Cauchy transforms (or good approximations) of 
/*/ and x. The Cauchy transforms of the elements x and f*f can be respectively computed as 
explained in equations (5) and (6). 

3.1. Numerical simulations. We compare the distributions obtained with our method and the 
empirical eigenvalue distributions of 20 realizations of 1000 x 1000 block-modified for the cases 
where x has a Wigner, Wishart and arcsine limiting distribution and the block transformation 


2) — ( 


f lla n T 15a22 — 25ai2 — 25a2i 36a2i 

36ai2 ll«n - 4a 2 2 


The numerical results are plotted in Figures 1,2, and 3. 



Figure 1. Eigenvalues of block-modified Wigner matrix (histogram) and asymp¬ 
totic distribution (solid line) computed with our method. 

Now we want to obtain explicit formulas. For this, we will some conditions on tp which allow 
us to express the asymptotic distribution of as the product of operators which are free with 
amalgamation over a commutative algebra. 

4. Freeness for block-modified random matrices 

We shall denote by {Ei\ r i=l a fixed basis of the complex Hilbert space C r . Also, we put Eij = 
E*Ej. To any linear map ip : M m { C) — > M n { C), one associates its Choi matrix [Cho75] 

m 

M n ( C) ( 8 > M m { C) 9 C = [<p ® id](H m ) = ^ <p{Eij) 0 E ij: (11) 

i,j= 1 

where the Bell (or maximally entangled) matrix is Bl m = oj m oj* n , with 

m 

oo m = ^E i ®E i eC m ® C m . 

2=1 


(12) 
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Figure 2. Eigenvalues of block-modified Wishart matrix (histogram) and asymp¬ 
totic distribution (solid line) computed with our method. 



Figure 3. Eigenvalues of block-modified randomly rotated arcsine matrix (his¬ 
togram) and asymptotic distribution (solid line) computed with our method. 


As in Section 2.4, we consider a sequence of self-adjoint unitarily invariant random matrices 
X d E M d { C) ® M m { C) that converges in distribution towards some random variable x E A having 
distribution (i. Define 

n m 

x% = [id ®<p](X)='£ T d ® E i:j )X d (I d ® E kl ) E M d ( C) ® M n ( C), (13) 

i,l=1 j,k= 1 

for some coefficients Cij k i E C, which are actually the entries of the Choi matrix of y? (see Lemma 
4.2) 


VI < i,/ < n, 1 < j,k < m, 


Cijki — {Eu ® Ejki C)' 
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In order to consider the general case where m ^ n, we embed the matrices Xj, as well as the 
matrix units Eij, Eki into the larger, square space M^(C) <8> M m+n { C), as follows: 






x d 

0 mxn Id 

0 nxm Id 

0 nxn ® Id 

0 mxm ® Id 

0 mxn ® Id 

0 nxm ® Id 

K 


n m 

EE Cijkl(Em+i,j Id) ' Xd ' ( Ek,m+l ® Id)- 

i,l =1 j,k =1 


Following Remark 2.4 (3), we have a nice description of the asymptotic distribution: 


Proposition 4.1. The sequence Xf } converges in distribution, as d —> oo, to an element xP E A 
of the form 

n m 

IP — ^ ^ ^ ('/] kl ni - i.j■!'('h’:m • l- 

i,l =1 j,k= 1 

where x and the family of matrix units {%} 7 y=i are free with amalgamation over the algebra 
( Pm,Pn ) C A. 

Here p m ,Pn are orthogonal the projections 

m m+n 

Pm, = ^ ^ £a and p n — ^ ^ &a , 
i =1 z=m+l 

with traces r(p m ) = m/(m + n), r{p n ) = n/(m + n). 

The distribution of x ism/(m+n)p+n/(m+n)5o. In particular, the random matrix X// converges 
in distribution, as d -A oo, to toe probability measure pP, which is such that the distribution of xP 
is m /(m + n)<5o + n/{m + n)p v . 


First, we exhibit the relation between the coefficients Cijki and the Choi matrix C of the linear 
map ip. 

Lemma 4.2. Consider the spectral decomposition of the Choi matrix (11) 

mn r 

c = ^2 X s v s v* = ^2 PtPu 

S=1 t= 1 

where the eigenvalues X s E M, are listed with multiplicity, v s E C n (g>C m , r is the rank of C, and 
Pt,t < r is the projector to the eigenspace with non-zero eigenvalue pt. Then, for all i,j,k,l = 
1 

r 

C-ijkl — {P%1 ® Pjk 5 C) — N ( X s (Ei ® Ej , V s ) (El ® Efc, v s ). 

S=1 

Proof. Use equation (13) to find c^ki = (Eu, ip(Ejk)). Since C = Ylij T{Pij)®Piji the first equality 
in the conclusion follows. The second one is obtained by using the spectral decomposition of the 
selfadjoint matrix C . □ 


Proposition 4.3. The block-modified random variable xP has the following expression in terms of 
the eigenvalues and of the eigenvectors of the Choi matrix C: 

x* = f*x®Cf, (14) 


/ = ^2 w * s ® v s e -4 <g> M nm (C ) 

S=1 


where 
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and the random variables w s G A are defined by 


n m 

W s = 'P, (Ei®Ej,v s)£m+i,j • 

i=l j =1 


(15) 


Proof. The proof follows directly from Lemma 4.2. For the left hand side of (14), we have 


i,l =1 j,k =1 


EE CijklCm+ijXCrn+l^ 

i,l =1 j',A;=l 


E EE ^s(Pi ® Pj i Vs) {Bl ® Ejci Vs)Grn-\-i t j3'€rn-\-l,k 

i,l =1 j,As=l s=l 


E> 

S=1 


s w s xw , 


Using the orthogonality of the eigenvectors u s , one can see easily that the right hand side of equation 
(14) gives the same expression as above, finishing the proof. □ 

For each element pt of the spectrum of C, let fi C {l,...nm} be the set of eigenvectors re¬ 
appearing in the eigenprojector P t : 

p t = J2 v i v *r 

jeJ t 

We also define Qt ■= YljeJt w j w j ^ 

We are going to compute first the distribution of the random variable 

y = (x®C)-ff *. 

The relation with the variable x v , in which we are ultimately interested in, is given by a normal¬ 
ization factor: 

\/k > 1, T [(, x<p ) k ] = nm ■ E[y fc ], 

In particular, the same relation holds for the \k-transforms: 

(16) 

Definition 4.4. We say that the eigenspaces of C are tracially well behaved if, for alli ±,..., < r, 

T(wj 1 w* 2 Q il ... Qi k ) = S jlj2 T{w jl w* jl Q il ... Q ik ) = S jlj2 T(w j ' i w*, i Q h ... Q ik ), (17) 
for every j \, j\, j? < mn such that ji,j[ G Ji for some i < r. 

Theorem 4.5. Consider a linear map ip : M m { C) —> M n { C) having a Choi matrix C G M nm (C) C 
A <%) M nrn (C) which has tracially well behaved eigenspaces (see Definition 4-4)- Then, the random 
variables x <S> C and ff* are free with amalgamation over the (commutative) unital algebra B = 
{Pm,Pn } <8> {C) generated by the projections p m ,p n from Proposition 4-1 and the matrix C. 

Proof. Clearly B is generated by p m ,p n and by the spectral projections P\,...,P r G M nm ( C) C 
A<8> M nm (C). Since x and (eij)ij are free with amalgamation over ( p m ,Pn ) (see Proposition 4.1), it 
is well known that (x) <g> M nm { C), (( eijfij) ® M nm ( C) C A <g> M nm ( C) are ( p m ,p n ) <S> M nm ( C)-free, 
with respect to the conditional expectation E := -0 <g) id, where ip (a) = T{p m a)p m + T(p n a)p n . Since 
ff* G {(eij)ij) < 8 > M nm ( C), x ® C and ff* are ( Pm,Pn ) < 8 > M nm (C)-free. By Proposition 2.8 and 
Corollary 2.9, to show £>-freeness, it is enough to show that the ( p m ,Pn ) < 8 > M nm (C)-moments of 
x <S> C and f f* restricted to B are also valued in B. 


= r ^ sc 
w x 1 ? 1 ^ y • 
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The case of x <g> C is immediate: for any i\, ..., ik-i G {m, n} and j 1 ,.. • , jk-i £ {1, ..., r}, we 
have 


E [(£ <g) C)(pi L <g> P h ){x <S> C) ■ ■ ■ (x <S> C)(pi k _ 1 <g> Pj k _!)(x < 8 > C)\ 

= ip^xp^x ■ ■ ■ xpi k _ x x) <g> CP h C ■ ■ ■ CPj k _ 1 C. 

Since x = p m xp m and p m p n = p n Pm. = 0, the first factor of the tensor product is zero unless 
i\ = ■ ■ ■ = ik~i = m, in which situation it is equal to r(x k )p m . It is also obvious that the second 
factor is an element of (C), so the whole expression is B- valued. 

Let us now treat the case of ff*. First, note that 

r 

ff* = Y w* h w h ' ® V h v* h , 
h,h '=1 

and, for any h, h', i, 


V* h ,PiV h = bh,h'1-hGJi 

Wh'PiW* h = 1 i=m.Wh'W* h . 


We develop 


E [//* (Pn ® P*) • • • 

(Pifc-1 ®-Pifc-i)//*] 

= E E 

K^h'fpi 1 ■ ■ ■ Pi k -l W h k W h' k ® • • • Pj k -!Vh k V* h ' k 

hi,...,h k = l 


K,..P k = i 



r 

= Y T ( W h 1 Qh ■ ■ ■ Qi k ~i w h' k )Pn ® v hiV* h ’ k 

hi,h' k =l 


r 

— Pn ® ^ ^ Cj /( G B , 
3 = 1 


(18) 

(19) 


( 20 ) 

( 21 ) 


for some constants Cj and the assertion follows. In the last equality above, we have used the 
hypothesis satisfied by the Pt (and thus by the Qt, see Definition 4.4) and the tracial property of 

T. ^ □ 


We would now like to investigate classes of Choi matrices C which satisfy the condition in 
Definition 4.4. 

We start with the following lemma: 

Lemma 4.6. For any s, £ G {1,..., r}, we have t(w s w*) = (rn + n)~ 1 5 s t. 


Proof. Using equation (15), we have 

n m 

t(w s w*) = E E i ® Eji •> Vs) (-®Z2 ® E j2 , Vt)T(Cm-\-ii, jl e j2,m+i2 ) 

il,i 2 =lji J2=l 
n m 

= Y Y ^*1 0 E h > v s) i E i 2 ® E h , vt)Si li2 Sj lj2 (m Pn)- 1 

il,i 2 =l jlj2=l 

= (m + nY l (v t i v s ) = (m + n) _1 <5^. 


□ 


Definition 4.7. T/ie Choi matrix C is said to satisfy the unitarity condition (UC) if, for all t, 
there is some real constant dt such that Qt = dt ■ 1. 
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Many Choi matrices associated to linear maps of practical importance satisfy the unitary condi¬ 
tion (UC), see, for example, Sections 7.3 and 7.4. 

Proposition 4.8. Any Choi matrix C satisfying the unitarity condition (UC) satisfies also the 
condition from Definition 4-4 ■ Moreover, dt = n~ 1 rkPt. 


Proof. Using (UC) and Lemma 4.6, we have 

T(w jl w* 2 Q il ... Q ik ) = d h d i2 ■ ■ ■ d ik T(w jl w* 2 ) = n~ l d h d i2 ■ ■ ■ d ik 5 jlj2 , 
which is enough to prove the statement in (17). For the second claim, note that 

dt = r{Qt) = ^2 T (. w j w j ) = n ~ 1 ^2 IKII 2 = n~ 1 Ti'P t . 
j&Jt jeJ t 


( 22 ) 


□ 


Remark 4.9. There are maps which satisfy the condition in Definition 4-4 but not the condition 
(UC). For example, the map ip : M n ( C) -A M m (C), given by: 

T ( A) — 'N ] atjEtjAEjt, 
i<n,j<m 

has a diagonal Choi matrix C = diag(an, a 21 ,..., a n i,ai 2 , ■ ■ ■, a nm )■ Hence the eigenvectors of C 
are just the standard basis elements ( Ej < 8 > F?fe)j,fe<n of C n . 

Therefore, the “matrization” Wjk of the eigenvector Ej ® E^ is simply the matrix unit e^. Thus, 
Qij = e u are projections and is does not satisfy (UC). However, we have that 

T ( w iiji w i2 j 2 Qi 3 j 3 ... Qi k j k ) = T(e hjl e j2i2 e i 3i3 ... ei k i k ) = n dj 1 j 2 ^j s i a+ i 5 (23) 

s<k 

which implies 17. 

By the discussion at the end of Section 2-4 we can actually compute the distribution of x v = 
Yli<m c a ® x a- Indeed, since the (xa) are free in An, the distribution /r of is just a convex 
combination of probability mesures /j, = I(ni + ■ ■ • + n n ), where each pLi is the distribution of a 
linear combination anxn + • • • + ai m x mm , where the (xjj) are free, identically distributed random 
variables, each having the same distribution as the free compression xn of x. 

We consider now the dual ip* of a linear map ip : M m ( C) -» M n ( C) acting on the blocks of a 
unitarily invariant random matrix. Recall that ip* : M n { C) —> M m ( C) is defined by the duality 
relation (with respect to the Euclidean, or Hilbert-Schmidt, scalar product) 

VX G M n (C),Y G M m (C), (X,<p(Y)) = {<T*(X),Y). 

The Choi matrices of ip and ip* are related in a simple way, see [FN15] and their eigenprojectors 
share similar properties. 


Lemma 4.10. Let ip : M m { C) —> M n ( C) be a linear map and p>* : M n ( C) -A M m ( C) be its dual 
with respect to the Euclidean scalar product. Then, 


Cp* — F n rn C^ F n m , 


where F U)m is the unitary operator acting on C n ®C m as F n)m x®y = y®x for iGC” and y G C m . 
In particular, if ip is hermiticity preserving, then C^ and Chave the same spectrum. 


Proof. The first statement is contained in [FN15, Lemma 4.3]. The second statement follows from 
the fact that the eigenvectors of C^* are obtained from the eigenvectors of C<p by coordinate-wise 
complex conjugation and swapping of the tensor factors. □ 
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5. Distribution of the product 

In the previous section, we have shown that freeness with amalgamation holds for our model, 
under some fairly general assumptions on the Choi matrix C. In this section, we use the freeness 
result to compute explicitly the distribution of the modified element we do however require 
some extra assumptions on the Choi matrix, in order to obtained a precise characterization of 
the modified distribution. If for the freeness result, the assumption that the eigenvectors C are 
well-behaved (see Definition 4.4) is enough, in the next result we require the stronger unitarity 
condition from Definition 4.7. 

Theorem 5.1. If the Choi matrix C satisfies the unitarity condition from Definition fil, then the 
distribution of x v has the following R-transform: 


S 

Rxr(z) = y diPiR x 
2—1 



(24) 


where pi are the distinct eigenvalues of C and ndi are ranks of the corresponding eigenprojectors. 
In other words, if p, resp. p v , are the respective distributions of x and x v , then 

p* = mu {D pi/n pf nd P (25) 

Proof. We start by computing the 0-valued ’I'-transforms of the elements x (g> C and ff*. For 
simplicity, we shall see the algebra B as a commutative subalgebra of M nm ( C) by identifying the 
elements 1 ® B and B, for any B E M nm ( C). We have 


= y E [(x <g> C • 1 <8> B) k 

k> 1 


yr(x k )(CB) k = * x (CB), 

k> 1 


where 4/ a .(-) is extended on B = (C) by (polynomial) functional calculus. The inverse 'F transform 
and the ^-transform are easily deduced from the above expression: 

= C'- 1 *<" 1> (£), S x ®c(B) = C- l S x {B) 

where (7 -1 is the pseudo-inverse of C , 

c - 1 = £ Pi l Pt. 


Let us now move on to the case of ff*. For a general element B = Y2t=i btPt, we have 


*//•(£) = 


fc>i 


ff* • 


i=i 


E E • • • ^-! E [//*(! ® Pti )//* • • • ff* (1 ® ■ B. 

k> 1 ti,...,tk-i = l 
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Using the unitarity condition (UC) and the relation (22), we get 


*//*(£) = E E • • • Vi E [//*(! ® ■■■//*(! ® • B (26) 

fc>l 


E E b ti--- b t k -i [J2 dt ^" dt k-i vhv *h n ls 


fc>i 


E E btdt 

fc>i \t=i 

n~ l B 


\h= 1 


fc-i' 


n 1 -B 


S 


(27) 

(28) 
(29) 


1 — n lr Tr B n —TVS’ 

where we have used the fact that dt = n _ 1 TVSj, see Proposition 4.8 (the above relation holds for 
small enough B ). The inverse T-transform and the 5-transform follow easily from the relations 
above: 




<-!>/ 


nB 


//* (^“i + Tr S’ 


= n 


B + I 


1 + TV B 

Since the elements x < 8 > C and //* are free with amalgamation over a commutative algebra, the 
5-transform of the product y and the inverse ^-transform are given by 

S y (B) = S x ® c (B)S ff *(B ) = nC7" 1 5 x (S) r |±E 

i s x(B) b 

y 


defats) = nC~ 


(30) 


1 + TVS 

Since we are only interested in the scalar distribution of the random variable y, we only need to 
invert the function on scalar elements, i.e. we need to solve the equation \EU _ 1 > (S) = al, 

for some real value a, and then we obtain T ty C (a ) = tr(S). In other words, using (16), we have 
^> x %(a) = Tr(S). Let B = Xa=i bj,Pi be a solution of the equation T^ 1 > (S) = al. On the 
eigenprojector Pi, this equation reads 


n biS s x c (bi 


= a. 


Pi 1 + TrS 

Let now Zi be such 6j = ZiR x c (zi). Using the formula zR(z) + S(zR(z)) = z and the previous 
equation, we obtain, for all i, 

Zi = —a( 1 + TrS). 
n 

Plugging this value into (a) = Tr(S), we get 

s s s 

T 8 x %(a) = TV(S) = ndA = E ndiZiRg(zi) = E d iPi a{l + Tr B)Rf \-a(l + Tr B) . 

i =1 2=1 2=1 

Replacing a by 1/(S*$, (a) + 1/a) and using TV(S) = T^(a) = aR s xV {a), we obtain a(l + TVS) = a. 
Thus, we finally get 


aR x %(a) = di P iaR s x 


2=1 


Pi ' 
—a 
L n . 


which is our conclusion. 


□ 


Remark 5.2. From the previous theorem, it follows that whenever the measure p is supported 
on the positive half-line, i.e p([0,oo)) = 1, and the Choi matrix C is positive semi-definite, the 
modified measure pP is also supported on the positive reals. This is a consequence of a linear 
algebra fact: C being positive semi-definite, the map p is completely positive, and thus, for any 
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positive semi-definite matrix X, the block-modified matrix X^ = [id® p\(X) is also positive semi- 
definite. However, neither converse holds: there exist matrices C with negative eigenvalues such 
that p v is not positively supported, although p is (see Section 7.3), and there exist measures p which 
are not positively supported such that p^ is positively supported (see the recent work [CHN15]J. 

6. Examples: specific distributions 

6.1. GUE matrices. Let us consider in this subsection the case of GUE matrices converging to 
a semicircular measure p sc of mean a and variance a 2 (see Section 2.4 for the definitions). This 
measure is characterized by the fact that only its first two free cumulants are non-zero: 

Rusc(z) = a + a 2 z. (31) 

Proposition 6.1. Consider a sequence of GUE matrices X d G M d [C) <S> M m (C) converging in 
distribution to the semicircular measure psc of mean a G M and variance a 2 > 0. Then, for any 
linear map p : M m (C) — > M n ( C) satisfying the unitarity condition (UC) from Definition f.7, the 
sequence of modified matrices Xv = [id <8> p\(X) G M d (C) <8> M n ( C) converge in distribution to 
another semicircular measure of mean aP and variance ( a^) 2 , where 


n 

a* = -JTr(C 2 ). 
n 

Proof. The proof follows from Theorem 5.1 and the form of the i?-transform for semicircular random 
variables (31). □ 

6.2. Compound Wishart matrices. Compound free Poisson distributions were introduced in 
[Spe98] as a generalization of the free Poisson (or Marcenko-Pastur) distribution. Given a positive 
measure p on the real line, the compound free Poisson distribution of parameter p is the unique 
(EB-infinitely divisible) probability measure on the real line having the moments of p as its free 
cumulants, see Section 2.4 and [HP00, Proposition 3.3.11]. 

As in the case of the usual free Poisson distribution, there is a very natural random matrix 
model for the compound free Poisson class of measures. Let A be the mass of p and consider 
the normalized probability measure p(-) = p(-)/\. Let X d G M d xs d ( C) be a sequence of random 
Ginibre matrices (i.e. random matrices with i.i.d. standard, complex Gaussian entries) such that 
Sd ~ A d. Consider also a sequence of selfadjoint matrices Dd G M Sd (C), independent from Xd (or 
deterministic) that converge in distribution to p. Then, the sequence of random matrices 

-X d D d X* d 

Sd 

converges in distribution to the compound free Poisson distribution 7r^; see [HP00, Proposition 
4.4.11], 

In the next proposition, we show that the class of compound free Poisson distributions is stable 
under linear block modifications. 

Proposition 6.2. Consider a sequence of compound Wishart matrices X d G M d ( C) ® M m ( C) 
converging in distribution to a compound free Poisson measure i r M , where p is a positive measure 
supported on the real line (not necessarily of unit mass). Then, for any linear map p : M m ( C) — > 
M n (C) satisfying the unitarity condition (UC) from Definition f. 7, the sequence of modified matrices 
X v = [id <g> p\(X) G M d ( C) ® M n ( C) converge in distribution to another compound free Poisson 
measure where 

S 

pV> = ^ndiD p ./ n [p]. 

1=1 
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In particular, if the limiting eigenvalue distribution of the random matrices Xd is a (usual) free 
Poisson measure of parameter X >0 (i.e. p = X8\), then the limiting eigenvalue distribution of the 
modified matrices X f is a compound free Poisson distribution of parameter 

S 

= A ^2ndid Pi / n , 

1=1 

which is Xmn times the empirical eigenvalue distribution of the rescaled Choi matrix C/n. 


Proof. Using Theorem 5.1 and the fact that the free cumulants of a compound free Poisson measure 
7 are the moments of p, we have (m p denotes the p-th moment of a measure) 


R x r 0) 


s 

} ^ dfpiRx 

i= 1 



s 

^ ^ dipiRix^ 
i=1 


'Pi 
. n 


z 


s oo 

^2dipi^2m p+ i(p) 

1=1 p =0 

s OO 

Y diPi — m P +1 ( D Pi/n M ) Z P 

i =1 pi p=0 

OO 

Y m P+i 

p =o 



□ 


Remark 6.3. The above result generalizes [BN12b, Theorem 3.1], in the case where the assumptions 
of both statements are satisfied. Note however that the unitarity condition (UC) from Definition 
4-7 is, in practice, much easier to verify that the set of conditions from [BN12b, Theorem 3.1], 
which are indexed by an integer parameter p E N. 


7. Examples: specific linear maps 

7.1. Unitary rotations of arbitrary random matrices. To warm up, let us start with the 
unitary rotation map p : M n (C ) —>• M n (C) given by <p(X) = UXU*, for some fixed unitary 
rotation U E U n . In this case, it is clear that for any (random) bi-partite matrix C) <S> M n ( C), 
one has 

= [id®p](X) = {I ®U)X(I ®U)*, 
so that and X have the same spectrum. 

Applying the formalism of our work to this case, we find that the Choi matrix corresponding to 
map <p has unit rank: C = uu* , where u E C n (g> C n is the “vectorization” of the square matrix U : 

n n 

U = Yj{ e iiU e j) e i ®Zj= Yh Uij e i ® e j- 

i,j=1 i,j =1 

Since ||u|| = y/n, the Choi matrix C has a unique non-zero eigenvalue, equal to n, of multiplicity 
one. Notice that C satisfies the unitarity condition (UC) from Definition 4.7, since this property is 
equivalent to the unitarity of the matrix U. Considering a sequence of matrices X^, converging in 
distribution to a limit p, we apply Theorem 5.1, verifying that pP = p. 
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7.2. The trace and its dual. In this subsection, we look at the following two maps: 

V’A : M m { C) -A C ^ Mi(C), VaPO = TV(T t X) 

J : M!(C) -A M n (C), J(z) = xJ„, 


where A E M m (C) is a self-adjoint matrix and i£l. Note that, for m = n, the maps J and Vh 
are dual to one other. The Choi matrices read, respectively, 

= Mi( C) ® M m (C) 

Cj = xl n ®1 £ M n (C) <g> Mi(C), 


and it is easy to check that they both satisfy the unitarity condition (UC) from Definition 4.7. 
In the case of the map J, the modified matrix is = zXd <8> I n , so its spectral distribution 
will converge to D x \y\, as predicted by Theorem 5.1. Using the same theorem, we can prove the 
following proposition. 

Proposition 7.1. Consider a sequence of unitarily invariant random matrices Xd E Md{ C) <8> 
M m ( C) converging in distribution to a probability measure y. Then, the sequence of matrices 
X^ A := [id <8> ip a] {Xd) E M^{ C) converges in distribution to the probability measure given by 

^ A = D Xl [y] ffl D X2 [/r] ffl • • • ffl D Xm [n ], (32) 

where Ai,..., \ m are the eigenvalues of A. In particular, when A = I rn , X^ 1 := [id < 8 > Tr](Xd) E 
Md( C) is the partial trace of Xd, and 

= ( 33 ) 


7.3. The partial transposition. Historically, the study linear modifications of block-random 
matrices was initiated by Aubrun in [Aubl 2 ], where he studied partial transpositions of large 
Wishart random matrices. The same random matrix model, in the asymptotical regime where one 
dimension is kept fixed, while the other grows, was analyzed in [BN12a], and later generalized in 
[BN12b]. 

Let us consider the transpose map, transp : M n { C) — > M n ( C), transp(X) = X T . Its Choi 
matrix is given by the flip operator F E M n (C) < 8 ) M n (C), which is the unitary matrix defined 
by F(x <E> y) = y <8> x for all x,y E C n . The spectral decomposition of the flip operator reads 
F = P + — P_, where P+, resp. P_, are the orthogonal projections on the symmetric, resp. anti¬ 
symmetric subspaces of the tensor product C n <S> C n . Note that the transposition map satisfies the 
unitarity condition (UC) from Definition 4.7; this is a consequence of the relations 


P+ = 


I + F 
2 


and 


P_ = 


I-F 

2 


and of the fact that both the identity and the flip have partial traces which are proportional to the 
identity. 


Proposition 7.2. Consider a sequence of unitarily invariant random matrices Xd E Md(C) < 8 > 
M n ( C) converging in distribution to a probability measure y. Then, the sequence of partially trans¬ 
posed matrices X^ := [id(g>transp] (Xd) E Md(C)®M n (C) converges in distribution to the probability 
measure y r given by 


y T = D 


1/n 


^EBn(n+1)/2 g ^EBn(n—1)/2 


( 34 ) 


where B is denotes following operation: u\ B V2 := v\ B D_ 1^2. 
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Proof. Since the eigenvalues of the Choi matrix F of the transposition map are +1, — 1, with 
respective multiplicities nd± = n{n + l)/2, nd ,2 = n(n — l)/2, we have from Theorem 5.1 


/ = (D 1/n pf^ +1 V 2 ffl (D_ 1/n vf n{n+1)/2 
= D l/n /J sin< ~ n + i y 2 a (n_ 1/ u) EBn ( n+i v 2 

= D 1/n L ®n{n+l)/2 B /i ffl„(n- 1 )/2 


□ 


In the case where the measure /i is a free Poisson measure of parameter A > 0, we recover 
Theorem 4.1 from [BN12a]: 


Corollary 7.3. Let Wd £ Md( C) <8> M n (C) 6e a sequence of Wishart random matrices converging 
in distribution to a free Poisson distribution of parameter X, tt\. Then, the partially transposed 
random matrices = [id ® transp] (IT^) converge in distribution to a rescaling of a free difference 
of free Poisson distributions of respective parameters Xn(n + l)/2 and Xn(n — 1 )/ 2 .* 

TTa D\/ n [^"An(n+1)/2 ^ Tt\n(n—1)/ 2 ] • 

For random projections, we recover a result proved in [FN15, Proposition 6.5]: 


Corollary 7.4. Let Pd £ Md(C) (g> M n (C) 6e a sequence of random projection operators converging 
in distribution to a Bernoulli distribution bt = (1 — t)8o + t5\, for some fixed parameter t £ (0,1). 
Then, the partially transposed random matrices Pj = [id 0 transp] (Pd) converge in distribution to 
a rescaling of a free difference of free additive powers of Bernoulli distributions: 


^A D\/ n 


, Eln(n-)-l)/2 


B b 


n(n— 1)/2 
t 


7.4. The reduction map. In [JLN14, JLN15], the authors study the reduction map , following 
a motivation originating in quantum information theory. The reduction map is another example 
of positive, but not completely positive application between matrix algebra, and its study has 
originated with the introduction of the reduction criterion [CAG99]. It is of particularly simple 
form: for X £ M n ( C), 

X red := red (A) = I n • TY(X) - X. 


Proposition 7.5. Consider a sequence of random matrices Xd £ M^(C) <g> M n (C) converging in 
distribution to ta probability measure /i. Then, the sequence of partially reduced matrices X^ ed := 
[id (8) red] (Xd) £ Md( C) ®M n (C) converges in distribution to the probability measure pL red given by 


pT ed = D 


1/n 


ffln 2 — 1 rn r) .. 
/i LJ U n —lfl 


(35) 


Proof. The Choi matrix of the reduction map red reads C re d = I—F. Its eigenvalues are respectively 
1 and 1—n, with multiplicities n 2 —1, resp. 1. Moreover, the eigenvector associated to the eigenvalue 
1 — n is a multiple of the Bell state (12), having a partial trace equal to the identity operator on 
C n . Hence, the Choi matrix satisfies the unitarity condition (UC) and, using Theorem 5.1, we have 


yf ed = (Dy^f™ 2 - 1 EH (D {1 _ n)/n ri 


jn 

D\/ n 


Ein 2 —1 m 7-i 

fl LLI U\— n pj 

Ein 2 —1 n n 

h i/i 


□ 


In the case where the measure p, is a free Poisson measure of parameter A > 0, we recover [JLN14, 
Theorem 9.1]: 
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Corollary 7.6. Let Wd E M^(C) <8> M n (C) 6 e a sequence of Wishart random matrices converging 
in distribution to a free Poisson distribution of parameter A, ir\. Then, the partially reduced ran¬ 
dom matrices Wf ed = [id <S> red^VFd) converge in distribution to a rescaled compound free Poisson 
distribution Di/ n n llriX , where 

Pn ,A — \(ri l)(5l “I - A(5i_ n . 

7.5. Generalized partial transposition and reduction maps. We generalize in this subsection 
the results from Sections 7.3 and 7.4, by considering the map ip : M n ( C) —> M n ( C) dehned by 

<p(A) = a A + /3Tr(A)I n + 'yA J . (36) 

Its Choi matrix reads 

Cip = o.Ll n + f 3 I n 2 + 'yF. 

The eigenvalues and eigenprojectors of C v are as follows: 

Ai = not + (3 + 7 P\ = W l Ll n 

A2 = /3 + 7 p 2 = 1 + 2 F ~ n~ 1 Ll n 

I - F 

A3 = /3-7 p -i = 2 ' 

Since all the partial traces of the eigenprojectors Pi, 2,3 are multiples of the identity operator, the 
Choi matrix C^ satisfies the unitarity condition (UC), and the following results hold (we leave the 
proofs to the reader). 


Proposition 7.7. Consider a sequence of random matrices Xd E M^(C) <g> M n (C) converging in 
distribution to ta probability measure p. Then, for the application p from (36), the sequence of 
the block modified matrices Xj := [id (g) p\(Xd) E M^(C) <g) M n ( C) converges in distribution to the 
probability measure pP given by 


pP Pa+(f3-\-j)/n [p] ® P{p+p/n 


E0n(n+l)/2-l 

P 


® P(fi—p)In 


EBn(n—l)/2 
/i 


(37) 


Corollary 7.8. Let Wd E M^(C) ® M n (C) be a sequence of Wishart random matrices converging 
in distribution to a free Poisson distribution of parameter A, tx\. Then, the block modified random 
matrices Wj = [id <g> p>\(Wd) converge in distribution to a compound free Poisson distribution 
where 


p = A 


^a+(/3+7)/n T 


n(n + 1) 


A n(n — 1) 

— *-J ^ 2 


The results above generalize the ones from Sections 7.3 and 7.4, with the choice of parameters 
a = (3 = 0 , 7 = 1 , resp. a = — 1 , f3 = 1 , 7 = 0 . 


7.6. Mixtures of orthogonal conjugations. Inspired by Weyl-covariant quantum channels (see, 
e.g. [DFH06]), we consider maps ip : M n (C) — > M n ( C) of the form 

n 2 

^p{A) = Y J ^iU i AU*, (38) 

i =1 

where a* E M and Ui E U n are orthogonal unitary operators: 

Tr (UiUj) = nSij. 

The Weyl-covariant channels from [DFH06, Section 4] are particular cases of the situation above, 
where {a*} is a probability vector and the Uj are the Weyl operators 
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for 1 , 1/6 hd and U, V being respectively the shift and phase operators 


UE r = E 


cc+1 


and 


VE X = exp 


/ 2irix\ 


\ n 


E x . 


Maps <p as in (38) are interesting because of the special form of their Choi matrix. Indeed, the 
Choi matrix reads 

n 2 

Cp - ^ ( 77 (}, /' 11/ i l: r j , 

i =1 

where the tq are orthonormal vectors in C n ® C n given by 


1 



E 


S,t= 1 


Ui(s,t)E s ® E t . 


Hence, the eigenvalues of Cp are {nctj}, they are simple, and the partial traces of the corresponding 
eigenprojectors are all equal to n~ 1 I n . Thus, the Choi matrix Cp satisfies the unitarity condition 
(UC). We leave the proofs of the following results to the reader. 


Proposition 7.9. Consider a sequence of random matrices Xg E M^(C) ® M n (C) converging in 
distribution to ta probability measure p. Then, for the application tp from (38), the sequence of 
the block modified matrices Xj := [id® (f\(Xd) E My{ C) ® M n ( C) converges in distribution to the 
probability measure pT given by 

p* = mt 1 D ai [p]- ( 39 ) 

Corollary 7.10. Let Wy E Mg(C) ® M n (C) be a sequence of Wishart random matrices converging 
in distribution to a free Poisson distribution of parameter A, tt\. Then, the block modified random 
matrices W ( f = [id® p\(Wd) converge in distribution to a compound free Poisson distribution 
where 

n 2 

p = \ ^2s a .. 

i= 1 


In the case where the coefficients a/ are all equal to 1, we obtain the map p>{A) = Ti^H)/^, thus 
recovering the special case a = 7 = 0, /3 = 1 from Proposition 7.7. 
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